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ABSTRACT 

Correct interpretation of acoustic travel times measured by time-distance he- 
lioseismology is essential to get an accurate understanding of the solar properties 
that are inferred from them. It has long been observed that sunspots suppress 
p-mode amplitude, but its implications on travel times has not been fully investi- 
gated so far. It has been found in test measurements using a 'masking' procedure, 
in which the solar Doppler signal in a localized quiet region of the Sun is artifi- 
cially suppressed by a spatial function, and using numerical simulations that the 
amplitude modulations in combination with the phase-speed filtering may cause 
systematic shifts of acoustic travel times. To understand the properties of this 
procedure, we derive an analytical expression for the cross-covariance of a signal 
that has been modulated locally by a spatial function that has azimuthal sym- 
metry, and then filtered by a phase speed filter typically used in time-distance 
helioseismology. Comparing this expression to the Gabor wavelet fitting for- 
mula without this effect, we find that there is a shift in the travel times, that 
is introduced by the amplitude modulation. The analytical model presented in 
this paper can be useful also for interpretation of travel time measurements for 
non-uniform distribution of oscillation amplitude due to observational effects. 

Subject headings: Sun: helioseismology — Sun: oscillations — Sun: sunspot 



1. Introduction 



Time-distance helioseismology (jDuvall et al.lll993l ) is a local helioseismological technique 
that has been used to study meridional flows, flows and sound speed perturbations beneath 
sunspots (e.g.. lKosovichev &: Duvallll 19971 : iGiles et al. l ll997l : lzhao. Kosovichev. fc DuvalfeoOlh 
It measures the time for a wavepacket to travel between any two points on the solar sur- 
face, by computing a temporal cross-covariance between the Doppler time series at the two 
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points. The travel time is then inverted to infer various properties that are useful to map 
local structures of the sun. These results are interesting, as they complement the results ob- 
tained from global helioseismology based on inversion of normal mode frequencies. However, 
many aspects of time-distance helioseismology are still not fully understood. In particular, 
it has been observed that sunspots suppress the p-mode amplitude appreciably, but its con- 
sequences on travel times and the properties derived by inverting them has largely remained 
unexplored. 

The commonly used procedure of measuring the phase an d group travel times of acous- 
tic waves is based a Gabor-wavelet fitting formula derived by iKosovichev &: Duvalll (119971 ) 
for cross-covariance of randomly excited oscillation modes of the quiet Sun (represented 
by a spherically symmetric model). Initially, this procedure included only a broad-band 
frequency filtering of the solar oscillation signal, centered at the peak of acoustic power. 
Later, it was modified by adding a phase-speed filtering procedure in order to isolate the 
first-bounce signals (direct wave s without addition al reflections from the surface) and to 
improve the signa-to-noise ratio (jDuvall et al.lll997l ). The phase-speed filtering is important 
for the analysis of acoustic-wave packets traveling short travel distances, e.g. less than 1 
heliographic degree (~ 12 Mm). For larger distances the travel times can be measured with- 
out the phase-speed filtering, and it is found that the use of the phase-speed filtering does 
not significantly affects the measurement results. For shorter distances the influence of the 
phase-speed filtering may be significant, and has to be taken account. 

In this paper we study the effects of the phase-speed filtering and local spatial suppres- 
sion of Doppler velocity signals fr om a quiet patch on th e Sun, on the travel times. In this 
'masking' procedure suggested by iRajaguru et al.l (120061 ) to simulate this effect in sunspots 
the oscillation signal of a quiet Sun region is multiplied by an inverted modulation function 
of spatial coordinates with azimuthal symmetry. This function is called a mask, and is not 
a function of time. One should note that while this procedure models a spatial modulation 
of acoustic power, it does not represent the modulation observed in sunspots, where the 
amplitude of acoustic changes due do several physical factors, such as reduced excitation, 
absorption, chang ing wave speed, and spectral line formation observing conditions. Recently, 
Chou et al.l (120091 ) attempted to quantify these contributions using observational data. The 
variations of the oscillation power in sunsp ots have not been explained by theory or simula- 
tions (e.g. IParchevsky fc Kosovichevl 120071 ). Nevertheless, it is interesting to investigate the 
effects of the variations on the helioseismic travel times by using the simple 'masking' model 
bearing in mind that it may not represent the the real situation in sunspots. One advantage 
of the 'masking' model is that it allows relatively simple analytical investigation. 



To understand the results of amplitude suppression or enhancement we theoretically 
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model the effect of masking by computing the cross- covariance of a masked signal that has 
been filtered by an appropriately designed phase speed filter. Analytical expressions for the 
cross-covariance are derived in terms of the mask, the filter parameters, the properties of 
the signal and the dispersive nature of the solar medium, when the the mask is azimuthally 
symmetric. This analysis will be useful to study the effect of masking on travel time mea- 
surements and the properties that are inferred from inverting the travel times. This will be 
especially valuable in artificially mimicking how the sunspots influence the travel times of p 
modes, and also other instrumental effects that corrupt the observed signal. 



2. Effect of phase-speed filtering on cross-covariance and travel times 



Consider a Doppler signal from a quiet patch on the Sun. Travel time m aps for this re- 

gion are computed by fitting the observed cross-covariance by a Gabor wavelet (jKosovichev fc Duvall 
19971 ). Now this region is masked by a spatial function to induce an artificial suppression in 
amplitude. One could alternatively enhance the amplitude in a similar manner. Maps for 
both mean and difference travel times are computed, also by fitting a Gabor wavelet. Taking 
the difference between the quiet and masked travel time maps, one sees app reciable shifts in 



the tr avel times around the masked region. This is illustrated in the paper of lRajaguru et al. 
(120061 ). It is generally observed that the difference travel times show appreciable shifts com- 
pared to the mean travel times in the masked region. 

In time-distance helioseismology we deal with acoustic waves near the solar surface, 
that are observed by measuring the line of sight Doppler velocity signal on the solar surface 
that has both radial and horizontal components of displacement. Without loss of generality, 
the line of sight direction is taken along the X axis. The signal is a sum of normal modes, 

— # — * 

and at a position R = (R, 9, (j)) on the solar surface and time t is written as, d x (R,t) 

E 



n i m d® lm (R,t) (jChristensen-Dalsgaardl 120021 . e.g.). For each mode, 



d nlm( R , t) = d nlm( R ) eX P ~ ^nlm)} 

where, the spatial part is given by projecting along the X-axis. 



d n lm( R ) - a rdmCrd{ R ) { Y l 

where, 
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(3) 



where, 9 is the co-latitude, (f> is the longitude, R is the solar radius, i = y— 1, the mode am- 
plitude a n i m , the phase a n ; m , and ^(R) and ^(R) are the radial and horizontal components 
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of the eigenf unctions respectively evaluated at R, and are therefore numbers. The spherical 
coordinate system is specified by unit vector e r in the radial direction r, and two unit vectors 
eg and in the horizontal directions 9 and <fi respectively. Each normal mode is specified 
by a 3-tuple (I, m, n) of integer parameters, corresponding angular frequency uj n i m , the mode 
amplitude a n i m , the phase a n i m . The integer / denotes the degree and m the azimuthal order, 
— I < m < I, of the spherical harmonic Y[ n {9, 0), which is a function of the co-latitude 9 and 
longitude <fi. These describe the angular structure of the eigenfunctions. The third integer n 
of the 3-tuple (l,m,n) is called the radial order. For a spherically symmetric Sun, all modes 
with the same n and I have the same eigenfrequency u n i, regardless of the value of m. 

In time-distance helioseismology we measure the travel time of wave packets by forming 
a temporal cross- covariance between the oscillation signals at two locations separated by an 
angular distance on the solar surface. To model this we represent the solar oscillations on 
the solar surface as a linear superposition of normal modes, that are band-limited in angular 
frequency uj. This is achieved by replacing a n i m ^ r nl {R) in equation <^ by the Gaussian 
frequency function Gi(u), which models the amplitude of the solar modes, 

Gi(uj) = h exp (-(«, - uj ) 2 /5lu 2 ) (4) 



where, bi is a coefficient of /, which is discussed in (INigam et al.ll2007l ). This function groups 
modes within a certain range of frequencies, which is controlled by the width 5u, about a 
central frequency u in the uo — I diagram. 

A phase speed filter is applied, and modes are selected from the uj — I diagram to 
construct the cross-covariance wave packet. It is specified by a Gaussian centered around a 
phase speed V p h and a width 8V p h as parameters, and is given by 

F P (V P ) = exp (-(V p - V ph ) 2 /5V p 2 h ) (5) 

where, the phase speed V p = u/L, L = + 1) = k h R, k h is the horizontal wave number. 
The role of the phase speed filter is to select waves with a small range of phase speeds, 
the range is specified by the width 5V p h- All these waves travel approximately the same 
horizontal distance on the solar surface, and sample the same vertical depth in the solar 
interior. Hence, it is crucial to select 5V p h appropriately so as to make the cross-covariance 
more robust, and hence be able to resolve the sub-surface structures in the Sun. 

Due to the band-limited nature of the oscillation amplitudes, only values of L which are 
close t o L Q = u) /V p contribute to the sum of equation (jTJ, and hence, following lKosovichev fc Duvall 



( 119971 ). we Taylor expand L about the central frequency lu , up to the first order: 



L(uj) = L[u + {to — 0J o )\ L(uj ) + - w ) (6) 
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The equation ([6]) can be written in terms of the group velocity U g = du/dL and phase 
velocity V p = lo/L, evaluated at u — u , and using the fact L{u ) = oo /V p , we have 

L(u) « £ + (1 - -L) Wo (7) 



Ug \V P U g 



Likewise, the phase velocity V p (L,uj) can be expanded about the point (L ,uj ) in the u — I 
diagram to yield 



2 



v2 



F l M.«,<)K<^\-Vj l \L-—) /5 f 

where, 5f = w 5V p h/V p h, and the filter width 5V p h is evaluated at (L ,lj ), and is a constant. 

The Taylor expansion is valid when the second order effects are small. These may not 
be small for small distances A, when the waves spend most of the time in the outer layers of 
the Sun. In these layers all the solar properties change abruptly and there are large gradients 
present, so higher order terms in the Taylor expansion should be retained. This could make 
the analytical calculation formidable. 

The cross-covariance ipj (Rx, R 2 , r) for the phase speed filtered Doppler signal 
d^R, 9, 0, t) = 1 m Fp(L, u> n i)d*i m {R, t) as a function of the time lag r is defined as 

T 

i}f p {R x ,R\,T) = iy dl(R u t)dl(R 2 ,t + r)dt (9) 


and involves integrating the product of the projected line of sight filtered Doppler signals 
at the two locations Ri = (R,9i,<f)i) and R2 = (R,9 2 ,(p2) on the solar surface over a time 
interval T, that is related to the period of the time series being cross correlated. Here we 
have replaced a n ; m ^ ; (i?) by the Gaussian frequency envelope function Gi{u) in equation (j4]). 

The cross-covariance from equation (Q is therefore, 

^(Ru^r) « Yl ^=^/f ^ Fp2(L '^ )G ' Ki)[/+KiT) + f - {UnlT)] (10) 

where, f + {ur) = cos(ur — LA + f — C) corresponds to the positive time lag and J-(ujt) = 
cos (ujt + LA — f + Q, corresponds to the negative time lag. Since cos is an even function, 
/_(c<jr) = /+(— out). This approximation is valid for large /, small A, such that LA is large. 
The phase term ( and the factor C\ are due to the horizontal component of the displacement , 



and depend on the location of the points R\ and R 2 , and are discussed in (INigam et al.ll2007l ). 
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The inner sum over u n i is replaced by an an integral over u, and we drop the negative 
lag term by extending uj to negative values to get 



/ 2 \ 

i>fl(R 1 ,R 2 ,T,V p )= J ^expl -2VJ(l-^ /5) j exp 



5to 2 



[UJ - L0 o f 



7T S 



cos(u;t— LA— C+ — 
(11) 



Evaluating the integral (IGradshteyn fc Ryzhiki Il994h . we get 



i>f p (Rl, R2, T, V p ) = A fp (5u, 8f, T, T g , T p ) COS {u) fp {j - T fp J + ^| (12) 

The shifted phase travel time 7/ h due to the phase speed filter and horizontal component is 
therefore, 



ph 



(Rg ~ Rp) T ph + 



l-Rn 



and the shifted frequency, Uf p = u (l — R gpe ). The amplitude scaling term is 



(13) 



.4 



5u e 



exp 



6u 2 e 



2J2 



8(m + s 2 



5oo 4 e 2 



(14) 



Summing equation (1T21) over phase velocities we get the final cross-covariance. 

iPf p (R u R 2 , t)=J2 V p) 
^ L V 7r A ' 



(15) 



The dimensionless quantities R g = r g — T ph /T ph and R p = t p — T ph /r ph represent the relative 
deviation of the group travel time r g = A/U g and phase travel time r p = A/V p respectively 
from the filter phase travel time, t p h = A/V p h- The filter widths appear in a dimensionless 
parameter e 2 = 5 2 /5uj 2 . The dispersive characteristics of the solar medium and the filter 
properties are related through the dimensionless parameter R gp£ = R g R p /(R 2 + e 2 ). 

Equations ( fT2"lTHj) provide a generalization of the Gabor wavelet fitting formula of 



Kosovichev fc Duvalll (119971 ) for the phase-speed filtering procedure. Obviously, the filter 



width 5V p h is very large, so that parameter e — > oo these equations are reduced to the 
standard fitting formula. 
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3. Effect of amplitude modulation on time-distance helioseismology 

measurements 

3.1. Cross-covariance for solar oscillations with spatial modulation, and 

travel-time shifts 

In this section, we derive a formula for the time-distance cross-covariance function in 
the presence of a localized amplitude masking. This provides estimates of the masking 
effect on the time-distance helioseismology measurements for various wave properties and 
observational parameters, including phase-speed filtering which is the major factor affecting 
the measurements. We consider a modulation function q(9) with azimuthal symmetry so as 
to simplify the analytical derivation. It can therefore be expanded as 



q(9) = ^^(cosfl) w J>V ^L6 COS { Ld ~ l) ( 16 ) 
i l 

where, the approximation is valid for large I, small 9, such that L9 is large (|Jacksonlll999l ). 
It should be noted that due to the azimuthal symmetry the mask function is independent of 
<j), and hence on m in this expansion in equation (fTEj) . Using orthonormality of Pi(cos8), we 
can compute the coefficient qi as 

qi = J q(9)P l (cos9)sm9 d9 (17) 

o 

Since the masking is carried out in a localized region we define the mask function that we 
apply to the signal as 



Q{9) = 1 + 8 q(9) = J2QiPi(^s9) « £qJ_!-oos (l9 - £) (18) 

i i 

where, s scales the mask function q(9), and is positive for enhancing the signal and negative 
for suppressing it. The coefficient Qi is 

7T 

Q l = ^i±l f Q(9)P l (cos 9) sin^ d9 (19) 



Masking takes place in the spatial domain and is independent of time. It consists of multi- 
plying the signal d x (R, 9, <fr, t) with the mask function Q(9). 
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The mask function Q(6) given in equation (fT8l) . and it is used to spatially modulate the 
signal d x (R,t) resulting in the masked signal dn(R,t), 

d x Q (R,t) = Q(6)d x (R,t) = [l+sq(6)}d x (R,t) = d x (R,t)+s q(6)d x (R,t) = d x (R,t)+s d x (R,t) 

(20) 

Where, d x (R,t) = q(9)d x (R,t) is the spatially modulated signal. The effect of masking is 
seen when we phase speed filter the masked signal. In the absence of a phase speed filter we 
do not observe any masking effect in the time-distance cross-covariance. 

Phase speed filtering of the masked signal leads to 

d% Q (R,t) = G l (u nl )F p (L,u nl )Q(6)d x nlm (R,t) (21) 

n,l,m 

The cross-covariance of the masked filtered signal is 

T 

r) = dl Q {R u t)dl Q {R 2 ,t + r)dt (22) 
o 

Substituting the expression for d x j Q (R,t) from equation (1211) into equation (122]) and using 
equation (1181) leads to 

^(R,, R 2 , r) = ^f p (R u R 2 , r) + s^f pJ jR u R 2 , r) + s ^f pqJp (Ru 4, r) + s 2 ^fjR u R 2 , r) 

(23) 

This is the final expression for the time-distance cross-covariance function with amplitude 
masking. It contains terms that result from the interaction of the phase speed filter and the 
mask function. The term ipj*(Ri, R 2 , r) is due to the phase speed filter alone, and does not 
contain the effect of the mask. The term ibf f is obtained by cross correlating the filtered 
signal with the modulated filtered signal, ipf pq j p is the cross covariance of the modulated 
filtered signal with the filtered signal. The last term ip* is got by cross correlating the 
modulated filtered signal at both the points R\ and R 2 . 

The cross-covariance from equation is therefore, 

T 

^%{R u R 2 ,t) = ^J dl{R u t)dl{R 2 ,t + r)dt (24) 
o 

^(R^R^r) = ^F p 2 (L,^)G2(^)cosK i r)< i (i? 1 ,il 2 ) (25) 
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Where, d^(i?x, R 2) = y-^rx co s(LA + Q is the m-averaged part of the spatial signal in the 
cross covariance (iNigam et al.l 120071 ). We have on combining the cosine terms, 



tZ&Mr) = £ ^T^^yf T, F p( L ^^M[f + (^nir) + Uu nl r)} (26) 



The phase factor ( is due to the horizontal component of the displacement, and depends on 
the location of the points being cross correlated. 

In equation (1261) we replace the dummy variable to n i by lu, and drop the negative time 
lag term by extending the sum over uj to negative frequencies, we get 



7T 



lot - LA + - - C, 



The inner sum in equation ( 1271) is written as 

(R u R 2 , r, V p ) = ]T F?(L, w)G?(a;) cos (cut - LA + - - C 

Multiplying with the mask function the different cross covariance are, 



(27) 



(28) 



r fpJ jRuR2,r) = - I dl{R x ,t)d% q {R 2 ,t + T)dt 



(29) 
(30) 



where 6 2 is colatitude of R 2 . 

Substituting the expansion for q{6 2 ) = ^ qiPi(cos9 2 ) into equation fl30|) we get 

iif pJpq {R u R 2 , t) = J2 Qi F p(L, co n i)Gf{u n i) cos{u n iT)P l {cose 2 )d x nl {R 1 , R 2 ) 

n,l 



(31) 



Substituting the asymptotic expansion for Pi(cos6 2 ) w J cos(L6 2 — |) into equation 
( l3Tj) and combining the cosine terms we get 



27TC, 1 /L 



vA^A27r V 2 V 7rL# 2 



"® £ F P ( L ' ^n^GfiuJni) COs(L6 2 --)[f + (uJ n iT)+f-(uJ n lT)] 



(32) 
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Combining the cosine terms in equation fl32l) and letting f + d 2 {u)r) = cos(ut — LA 2 - — C) , 
/-«o(wr) = cos(^r+LA 2 _+C) = f+da(—WT), f+p2^r) = cos(W-LA 2+ -C + f ), /- p2 (^) = 
U p2 (-lut), A 2+ = A + A 2 _ = A - 2 , we get 

V£ /pg 0Ri, 4, r, V p ) = £ F p 2 (L, WnJ )Gf K0(/ +£22 + /- d2 + / +p2 + /- p2 )/2 (33) 

In equation (1331) the functions / are evaluated at c<j n z. The sum can be converted into an 
integral over uj as before after dropping the negative lag terms /_<j and f- p , and extending uj 
to take negative values we get a sum of two Gabor wavelets. This shifts, the various travel 
times to r p2+ = A 2+ /% = T p + r pq2 , T pq2 = 2 /V p , r g2 = A 2+ /U g = T g + r gq2 , r gq2 = 6 2 /U g and 
r p h 2 = A 2+ /V p h = r p h + r p h q2 , T p h q2 = 2 /V p h. These shifts in travel times are related to the 
mask position. Hence, R g , R p and R gpe change to R g2 +, R p2 + and R gpE 2+ respectively, with 
the usual definitions. Similarly we have travel times for A 2 _: r p2 _ = T p — r pq2 , r g2 _ = r g — T gq2 , 
T ph 2- = T ph - r phq2 . Therefore, 




2l l>fp,fp*(Ri> ^2, r, V p ) = A fp (r g2 _, r p2 J) cos {uj fp2 _(r - r fph2 _)} + (35) 

A fp (r g2+ , r p2+ ) cos {uj fp2+ (r - r fph2+ ) + f } 

Similarly for the other terms 

T 

1%,j,{Ri,ik,T) = dl q {R^t)dl{R 2 ,t + r)dt (36) 



< aj/p (i?i, R2, r) = J2 J&^v! Ql ^ Fp2(L ' ^GK^PiicoseJiUi^r) +/_Krr)] 

Vp V UJ„l 

(37) 

Simplifying in a similar manner the inner sum of equation (|37|) is, 

2 C/^> ^, r ' ^ = r ^-) cos ( r - V-)} + (38) 

%( r 9i+> 7- pl+ ) cos {^/ pl+ (r - r fphl+ ) + I} , 

where subscript 1 in the travel times corresponds to the mask position 9\ of Therefore, 
C/ f *> r) = £ ^^vS 9 ^ /p9 ' /p( ^' ^' T ' ^ (39) 
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Finally, 

T 

^fjR u R 2 ,r) = iy dl^R^dl^t + r)dt (40) 



(41) 

Combining the various cosine terms in the inner sum of equation ( 14TI) we obtain, 
^fjRuR^Vj = ^F p 2 (L, W )G z 2 M[/ +1 + /_ 1 + / +2 + /_ 2 + / + 3 + ^3 + / + 4 + /-4]/4, (42) 

where /+i(wr) = cos(u;r - LA+i 2 _ - f - C), /-iM = / + i(-wr), /+a(wT) = cos(^r - 
LA +12+ + f - C), /_ 2 (wr) = /+ 2 (-wr), Z+aH = cos(u;r - LA_ 12 _ + f - C), /-s(wr) = 
/ +3 (-^r), /+4(wr) = cos(^r-LA_ 12+ + f-C), /- 4 (^r) = /-h(-wt), A +12+ = A+(0!+0 2 ), 
A +12 _ = A - {6 l + 2 ), A_ 12 _ = A - (8i - 2 ), A_ 12+ = A + {9 l - 6 2 ). The sum can be 
converted into an integral as before, and the resulting expression is a sum of four Gabor 
wavelets given by 





gi(Ri, R 2 , t, V p ) = A fp (T g+12 -, r p+ i 2 _) 


cos {uj fp+12 _ 


(r- 


T/ph + 12- 


)} 


(43) 




g 2 (R u R 2 , t, V p ) = A fp (r g+12+ , t p+12+ ) 


COS {^/p+12+ 


(r~ 


T/ph+12+ 


)} 


(44) 




g 3 (Ri, R2, r, V p ) = A fp (r g _ 12 _, r p _ 12 _) 


cos {uj fp _ 12 _ 


(r~ 


T /ph-12- 


)} 


(45) 




g^Ri, R 2 , r, V p ) = A fp (r g ^ 12+ , r p _ 12+ ) 


COS {^/p-12+ 


(r~ 


' /pfc-124 


)} 


(46) 


This shifts, 


the various travel times to r +pl2 _ = 


A +12 _/y p 


= t p 


~~ T P<?1 


— ?pg2j r H 


-5,12- = 


A +12 _/C/ 9 


— T g ~ T gqi ~ T 9<& an< i T +ph\2- = A_ 


{■12- /Vph = 


Tph 






These 



shifts in travel times are related to the mask position. Hence, R g , R p and R gpe change 
to i? +9 i 2 _, R +p \ 2 - and R +gp£ i 2 - respectively, with the usual definitions. Similarly, the other 
combinations can be defined. 



Therefore, 



t/jf (Ri,R 2 ,t) = Y -^2= i=* 2 ^f (Ri,^2,r,V p ) (47) 



The cross covariance ijjf (Ri, R 2 , r, V p ) is the sum of the four Gabor wavelets that depend 
on the masking parameters and is given by 



W fpq {Ri, R 2 , r, V p ) = g 1 {R 1 , R 2 , r, V p ) +g 2 (Ri, R 2 , r, V p )+g 3 (R u R 2 , r, V p )+g 4 (R u R 2 , r, V l 



v) 



(48) 
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Similarly, the cross covariance for the masking function Q(6) can be found with q\ repalced 
by Qi, and is given by 



where the coefficient Qi can be calculated from equation (1191) for the mask function Q{9). 
These equations which are nothing but the sum of Gabor wavelets can be used for the fitting. 
Comparing the different equations we find that the form of the Gabor wavelet is retained 
when a mask function with azimuthal symmetry is used. However, masking introduces shifts 
in the group and phase travel time by modifying the angular distance A. 

The formula in equation (J3BD is the masked cross-covariance. There are few extra 
parameters that represent the shifts in the phase and group travel times, due to the masking 
process that depend on the angular positions of the mask function. The coefficients qi 
an d Qi represent the fu nctional form of the mask functions. Similar dependence is seen 



in (IRajaguru et al.ll2006l ). A detailed numerical investigation of the amplitude modulation 
effects is outside the scope of this paper. In Figure (U we just give an example of the expression 
in equation fj48|) plotted for a particular mask position and compared with equation (Ti~2l 
without masking. We observe shifts in the cross covariance and hence the travel times change 
due to the masking process. To investigate to effect of the shape of the mask function, 
different values of Qi need to be included when computing the sum in equation fj4"9~l) . In 
this derivation we assumed that the modulation function exhibits azimuthal symmetry. For 
a general modulation function, the analytical approach becomes difficult, and numerical 
methods must be employed. 

In the model presented here, we assumed that the solar p modes have a narrow Gaussian 
amplitude function, that is peaked at u . Using t his fact we evaluat e the phase shift factor 



due to the horizontal component at this frequency (INigam et al.ll2007l ). Hence to a first order 
approximation, the phase shift is not affected by the masking procedure. In order to see the 
effect of masking on the horizontal component, we have to retain the frequency dependence 
of the phase shift factor when evaluating the integral in u. This will make the analytical 
approach intractable and the integral will have to be evaluated numerically. 



3.2. Phase-speed filtering and amplitude modulation (masking) do not 

commute 



In the previous section the signal d x (R,t) was first masked by Q(9) and then phase 
speed filtered by F p (L,uj n i). The signal is, 



d% Q (R,t) = J2 GW^^QW^t) (50) 



Expanding Q{6) from equation (fT8l) leads to 



(51) 



n,l,m 



We immediately see that this corrupts the eigenfunction due to the Pi(cos6) term in the 



filtering are not interchangeable. Also we observe that phase speed filtering followed by 
masking results in the cross-covariance being just scaled in amplitude by the product of the 
mask functions Q(0i)Q(9 2 ), the travel times being unaffected by the process of masking. This 
is because masking changes the mode eigenf unctions due to the presence of an additional 
term Pi(cos6) in equation (jSTl) . Of course, the modified eigenfunctions no longer satisfy the 
original wave equation for solar oscillations. Thus, this procedure cannot be considered as a 
physical model of the amplitude variations observed in sunspots. Our calculations show that 
the phase-speed filtering procedure of the masked signal leads to spurious shifts in travel 
times due to mixing of different I modes by the phase-speed filter. These shifts can occur 
if the masking is attributed to instrumental effects. However, this does not mean that the 
amplitude modulation due to the physical effects in sunspots results in similar effects in the 
travel-time measurements. This must be studied using realistic MHD models of sunspots. 




T 



^(RiMr) = \\ d* Qp {R u t)d) Qp {R 2 ,t + r)dt (52) 
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3.3. Effect of Gaussian frequency filtering 

Masking is observed only when the data is phase speed filtered after multiplying the 
signal with the mask function. In this section we show that just using a Gaussian frequency 
filter without a phase speed filter leads to no masking. Consider a masked signal dq(R,t) 
that is filtered by a Gaussian frequency filter G{ui) = exp (—(a; — uj ) 2 / 5uj 2 ). The mask 
function Q(9) is independent of /, and since the Gaussian frequency filter G(u>) does not 
depend on /, unlike the phase speed filter, the mask function Q{9) can be pulled out of the 
sum in equation (1531 

d % Q A f ) = E G(u nl )Q(e)d* nlm (R,t) = Q(9)d%(R,t) (53) 

n,l,m 

— » — # 

where, d x j g {R,t) = X] n z m ^(^)^nZm(-^;^) is the signal filtered by a Gaussian frequency 
filter. We see that the frequency filtered signal <^/ fl (-R, t) is scaled by the mask function Q(9), 
hence the cross-covariance ipj x Q (Ri, R 2 , r) of dJ gQ (R,t) is 

^fjR 1 ,R 2 ,r)=Q(9 1 )Q(9 2 )^(R 1 ,R 2 ,r) (54) 

where, ipj x g (Ri, R 2 ,t) is the cross-covariance of the Gaussian filtered signal dJ g (R,t). We 
see that the effect of masking on the travel times is not observed in ipj x Q (Ri, R 2 , r). It only 
undergoes a scaling of its amplitude by Q{9i)Q{9 2 ). So the effect of masking is not observed 
when the phase speed filter is absent, and we apply only a Gaussian frequency filter to the 
data. 



4. Conclusion 

In this paper we give an explanation as to why the localized spatial suppression or en- 
hancement (masking) of the acoustic signal followed by phase speed filtering can appreciably 
shift the measured travel times, and cause systematic errors in time-distance inversions. Us- 
ing only a frequency filter or reversing the operation of suppression (or enhancement) and 
phase speed filtering does not shift the travel times, but merely scales the cross covariance. 
Hence, the operations of masking and phase speed filtering are non-commutative. To explain 
this we develop a model and derive a new analytical formula for the cross-covariance in terms 
of the masking parameters, when the mask function has azimuthal symmetry. The reason 
for this is due to the fact that masking changes the spatial mode eigenf unctions, which when 
phase speed filtered, leads to a mixing of spatial modes, and is responsible for the travel 
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time shifts. This may be useful to mimic amplitude changes in sunspots due to pure obser- 
vational reasons, such as caused by altering the height of formation of the spectral line used 
for helioseismology measurem ents, and also to c orrect for the shifted travel times due to such 
effects. However, contrary to iRajaguru et al.l (120061 ) suggestion, the procedure of masking 
oscillations of the quite Sun cannot model the amplitude reduction in sunspot due to physical 
mechanisms (e.g. changes in emissivity or wave absorption), and thus their recommenda- 
tions of correcting the amplitude reduction by reversed masking may cause artificial shifts 
in observed travel times, and thus must be taken with caution. 
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Fig. 1. — An example of the cross- covariance function normalized by its maximum value as a 
function of the time lag. The solid line represents no masking that is computed from equation 
(|T2|) . while the star (*) line represents the effect of masking, calculated from equation (|48|) . 
It can been seen that there is a shift in the travel times as a result of masking. 



